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Abstract 

This paper contains three parts relating to the nucleon spin structure in a sim¬ 
ple picture of the nucleon: (i) The polarized gluon distribution in the proton is 
dynamically predicted starting from a low scale by using a nonlinear QCD evo¬ 
lution equation-the DGLAP equation with the parton recombination corrections, 
where the nucleon is almost only consisted of valence quarks. We find that the 
contribution of the gluon polarization to the nucleon spin structure is much larger 
than the predictions of most other theories. This result suggests a significant orbital 
angular momentum of the gluons is required to balance the gluon spin momentum; 
(ii) The spin structure function (f { of the proton is studied, where the perturbative 
evolution of parton distributions and nonperturbative Vector Meson Dominance 
(VMD) model are used. We predict asymptotic behavior at small x from lower 
Q 2 to higher Q 2 . The results are compatible with the data including the HERA 
early estimations and COMPASS new results; (iii) The generalized Gerasimov-Drell- 
Hearn (GDH) sum rule is understood based on the polarized parton distributions of 
the proton with the higher twist contributions. A simple parameterized formula is 
proposed to clearly present the contributions of different components in the proton 
to T\(Q 2 ). The results suggest a possible extended objects with size 0.2 — 0.3 fm 
inside the proton. 

PACS number(s): 12.38.Cy, 12.38.Qk, 12.38.Lg, 12.40.Vv 
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1 Introduction 


A precise determination of the polarized gluon distribution Sg(x, Q 2 ) is important in 
order to understand the spin structure of the nucleon. However, the direct measurement 
of the polarized gluon distribution in the nucleon is difficult. In the global analysis of 
the polarized lepton-nucleon deep inelastic scattering (DIS), the distribution Sg(x,Q 2 ) 
is extracted from the spin structure function gi(x,Q 2 ) through scaling violation as a 
higher-order effect of quantum chromodynamics (QCD). Unfortunately, such indirect de¬ 
termination of 5g(x,Q 2 ) is affected by large uncertainties because of the limited range in 
momentum transfer at fixed Bjorken-x and almost entirely arbitrary input gluon distri¬ 
bution. In fact, the data suggest that such global fit with either positive or negative input 
gluon distributions provides equally good agreement. Recently, a high-precision measure¬ 
ment of the mid-rapidity polarized proton-proton (p-p) collisions stringently constraint 
the polarized parton distribution functions mentioned above. The analysis of NNPDF 
collaboration [1] found an evidence for possible larger gluon spin distribution, which is 
against the common belief that it is rather small. 

In Sec. 2 of this work we use a QCD dynamic model of the parton distributions 
to predict the polarized gluon distribution in the proton without unknown input gluon 
distribution. Our model imagines that all gluons in the nucleon are radiated from the 
intrinsic quarks beginning at a low resolution scale. Thus, we can predict the radiative 
(unpolarized and polarized) gluon distributions provided the initial quark distributions 
are fixed. Such quark model was early proposed by [2,3,4] in 1977, it was improved in 
our previous work [5,6], where the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
equation [7,8,9] with the parton recombination corrections is used to reproduce the un¬ 
polarized parton distributions of the nucleon at Q 2 > 1 GeV 2 [10,11,12], Since the similar 
corrections of the parton recombination to the polarized DGLAP equation have been 
proposed in work [13], we can use these two modified evolution equations to predict the 
polarized gluon distribution in the proton dynamically. 

Comparing with the global analysis via scaling violation, the polarized gluon distribu¬ 
tion in this work is determined directly by the observed spin structure function gi(x, Q 2 ). 
We find that the contribution of the gluon polarization to the nucleon spin structure is 
surprisedly large, which is in excess of the previous estimations in theory. The reasons are 
as following: (i) The shadowing effect of the gluon recombination in the evolution of the 
polarized parton distributions is weaker than that in the unpolarized case since gdg « g 2 
at small x. Therefore, much more strong polarized gluons are emitted by quarks inside 
the polarized proton through a long evolution length from g 2 to Q 2 > 1 GeV 2 ] (ii) The 
positive contribution of the polarized gluon recombination, which is opposite to that of 
the unpolarized gluon recombination enhance the accumulation of the gluon helicity. The 
QCD evolution of the parton distributions begins from a low bound state scale g 2 not 
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only dynamically determine the polarized gluon distribution, but also exposes a novel 
spin-orbital structure of the nucleon in the light-cone frame, where the nucleon spin crisis 
has a possible explanation. 

Concerning the spin structure function, recently COMPASS experiment at CERN 
collected a large number of events of polarized inelastic scattering off the protons with 
very small x [14], The preliminary analysis of these data combining with the previous 
experiments [15,16,17,18], showed non zero and positive asymptotic structure function 
g\. In these fixed target experiments the low values of x are almost reached by lowering 
the values of Q 2 . The knowledge of the nucleon spin structure function gi(x, Q 2 ) at low 
Q 2 and small x is particulary interesting, since it is not only an important information 
to resolve the ”proton spin crisis”, but also provides us with a good place to study the 
transition from the perturbative research to the nonperturbative description of the proton 
structure. 

In Sec. 3 we try to study the behavior of g\ at small x but in the full Q 2 range. As 
we know that the structure functions of the nucleon are mainly constructed of the parton 
distributions at Q 2 > lGeV 2 , while the nonperturbative contributions to the structure 
functions become un-negligible at Q 2 C 1 GeV 2 . A key question is what components 
construct the spin structure functions of the proton at such low Q 2 1 Particularly, do the 
parton distributions and their pQCD evolution still play a role or not? To answere these 
questions, we discuss the application of the DGLAP equation with the parton recombina¬ 
tion corrections at low Q 2 in detail. We point out that the isolation of the contributions of 
the vector meson is necessary for keeping the factorization schema of the polarized parton 
distributions at low Q 2 . We find two different asymptotic behaviors of g\ at x < 10~ 3 : 
nonperturbaive behavior ~ x^ 1 at Q 2 < 1 GeV 2 and perturbative drop at Q 2 > 3 GeV 2 . 
We predict the translation of g l { at small x from lower Q 2 to higher Q 2 . The results 
are compatible with the data including the early HERA estimations and COMPASS new 
results. We point out that the measurements at different x with different values of Q 2 
in the fixed target experiments mix the complicated asymptotic behavior of gY The pre¬ 
dicted strong Q 2 - and ^-dependence of g\ at 0.01 < Q 2 < 3 GeV 2 and x < 0.1 due to 
the mixture of nonperturbative vector meson interactions and the QCD evolution of the 
parton distributions can be checked on the next Electron-Ion Collider (EIC). 

There is particular interest in the first moment ri(Q 2 ) = clxgi(x, Q 2 ) of the spin 
structure functions gi(x, Q 2 ), which has been measured from high Q 2 down to ~ 0 GeV 2 . 
The goal to obtain universal expressions describing Ti((5 2 ) at any Q 2 is an attractive task 
for both theoretical and phenomenological point of view. In theory, Pi(0) is constrained 
by the Gerasimov-Drell-Hearn (GDH) sum rule [19,20]. In Sec. 4 we try to expose 
the partonic structure in the GDH sum rule. Since we have known the contributions 
of g PGLAP+ZRS an d gY MD , one can expose the properties of TY t (Q 2 ) after subtracting 
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these two contributions from the experimental data about rf(Q 2 ). This opens a window 
to visit higher twist effects at low Q 2 in the nucleon structure. We proposed a simple 
parameterized form of T±(Q 2 ). We find that the negative twist-4 effect dominates the 
suppression of rf(<5 2 ) at Q 2 < IGeV 2 , while both the twist-4 and twist-6 effects have a 
dramatic change of Tf(Q 2 ) at Q 2 ~ IGeV 2 , which suggest a possible extended objects 
with size 0.2 — 0.3 fm inside the proton. 

Finally, following the above mentioned discussions, we will give a summary in Sec.5. 
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2 Dynamical determination of gluon helicity 
distribution in the nucleon 

2.1 Nonlinear polarized QCD evolution equation 


We use f + (x,Q 2 ) and f-{x,Q 2 ) to refer to parton (/ = q,q,g ) densities with positive 
and negative helicity which carry a fraction x of the nucleon momentum. The difference 
8f{x, Q 2 ) = f+(x, Q 2 ) — f-(x, Q 2 ) measures how much the parton of flavor f remembers its 
parent’s nucleon polarization. The spin averaged parton densities are given by f(x, Q 2 ) = 
f+(x,Q 2 ) + f-(x,Q 2 ). 

The spin-dependent QCD evolution equation of parton distributions with parton re¬ 
combination corrections was first derived by Zhu, Shen and Ruan (ZRS) in [13], it reads 


Q 


oi s {Q 2 ) r 1 dyx 


2 dx8q v (x, Q 2 ) 


dQ 2 


f --y8q v (y,Q 2 )AP qq i-), 
Jx v v y 


2vr Jx y y 


for flavor non-singlet quarks; 
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dQ 2 
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~ 4? tR?q 2 L' ^ xAP 99^9( x ^y)[y9(y,Q 2 )y^g(y,Q 2 )\ 

+ 4tt fpq 2 // 2 ^ xAp ^9( a:: ’ ?/ )[^( ?/ ’ ( 5 2 ) ?/ ^( ?/ ’ ( 5 2 )]’ x - V 2 )> 

2 dx8g(x,Q 2 ) 

^ dQ 2 

^ f d i X v [y £ ^ ^ 2 ) Ap <) + g 2 ) ap S3 (^)] 

Z7r j x y y i=1 y y 

+ / f s l?J 2 / 7 —xAP gg ^ g (x,7/)[7/c/(7/,g 2 )7/(5c/(?/,g 2 )], (i/ 1/2 < x < 1) (2.1.3) 

for gluon, where the factor l/(47rP 2 ) is from the normalization of two-parton distribution, 
R is the correlation length of two initial partons, the linear terms are the standard spin- 
dependent DGLAP evolution and the recombination functions in the nonlinear terms 
are (see appendix A) 
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27 (2 y — a:)(—20 y 3 + 12 y 2 x — x s ) 


64 
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y o 
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48 
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(2.1.4) 


(2.1.5) 


The spin structure function g\ at leading order (LO) and Q 2 > lGeR 2 is written as 


9 i(x, Q 2 ) = \ Y, e i [dqi( x , Q 2 ) + Sq^x, Q 2 )], (2.1.6) 

i 

where e t is the electric charge of the (light) quark of flavor i, i — u, d, s. 

The solutions of Eqs (2.1.1-2.1.3) are coupled with the spin-averaged evolution equa¬ 
tions, which are 


for valence quarks; 
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2 dxq v (x, Q 2 ) 
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cj(Q 2 ) 

4itR 2 Q 2 
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(2.1.7) 
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f A^R 2 Q 2 L/2 ~ XjP 99^( x ’^)[^(^’ ( 5 2 )] 2 ’( i -/’ x - V 2 ), 


2 dxqj{x,Q 2 ) 

^ dQ 2 

= ^ f 1 d ^ x -[y^)P qq {-) + y 9 M*)P qg ( x -)] 

2n J x y y y y 

+ AnR 2 Q 2 L/2 ^f xP 99^i( x ^y)\.y9(y.Q 2 )] 2 Aif V 2 < ® < 1), (2.1.8) 

for sea quarks; 
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rfg 2 
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AnR 2 Q 2 Jx y 


<*s{Q) r d y 


47 tR 2 Q 2 Jx1 2 y 


f —xPgg^g{x,y)[yg{y,Q 2 )] 2 , (if x < 1 / 2 ), 
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n ~ 2 dxg(x,Q 2 ) 
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^ y y 2^1 y y 

+ AnR 2 Q 2 Jx /2 ^ xP ^»( X ’ S/ )[^( J/ ’ < 5 2 )] 2 »(*^ i/ 2 -®- !) ( 2 - L9 ) 

for gluon, where the linear terms are the standard DGLAP evolution [7,8,9] and the 
recombination functions in the nonlinear terms are [10,11,12] 
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2.2 Dynamically radiative polarized gluon distribution 

We focus on the gluon distribution, its evolution is dominated by the valence quark 
distributions. Therefore, the corrections of the asymmetry sea distributions in the nucleon 
are neglected in this work. 

As we know, there are many effective QCD theories which describe the nucleon as a 
bound state of three quarks in its rest frame. The distributions of these quarks in the light- 
cone configuration have a similar character as the valence quark distributions observed 
at high Q 2 in DIS. Besides, the QCD evolution equation shows that either the second 
moment (i.e., the average momentum fraction) of the unpolarized gluon distribution or 
the first moment (i.e., the total hclicity) of the polarized gluon distribution increase 

as Q 2 increasing. A natural suggestion is that all partons (valence quarks, sea quarks 

and gluons) at high Q 2 scale are evolved from three initial valence quarks via QCD 
dynamics. Such idea was first proposed in 1977 for unpolarized parton distributions by 
[2,3,4], They assumed that the nucleon consists of valence quarks at a low starting point 
li 2 ~ 0.064 GeV 2 (but is still in the perturbative region a s (ii 2 )/2n < 1 and /i > A qcd), 
and the gluons and sea quarks are produced at Q 2 > /x 2 using the DGLAP equation. 
However, such natural input is failed due to the too steep behavior of the predicted parton 
distributions at small x since a long evolution distance from /x 2 to Q 2 > 1 GeV 2 . Recently 
the above naive idea was realized in the unpolarized DGLAP equation with the parton 
recombination corrections at LO approximation in [5,6], where the input distributions at 
/x 2 = 0.064 GeV 2 were extracted through fitting Ef ’ n (x, Q 2 ) and they have been fixed as 

xu v (x, IX 2 ) = 2A.3x 1 - 98 (1 - x) 2m , (2.2.1) 

xd v (x, fi 2 ) = O.lOa; 1 ' 31 ^ — a;) 3 ' 8 , (2.2.2) 

while 

g(x, fi 2 ) = 0, q t (x, li 2 ) = q t (x, /x 2 ) = 0, (2.2.3) 

and the parameters in Eqs. (2.1.8) and (2.1.9) are Aqcd = 0.204 GeV and R = 4.24 GeV~ l . 
We plot these input distributions in Fig.l. 

Similarly, using the polarized DGLAP equation with the parton recombination cor¬ 
rections Eqs. (2.1.1)-(2.1.3) and combining Eqs. (2.1.7)-(2.1.9), we fit g\' n (x, Q 2 ) with 
the data [21] in Fig.2, and extract the input polarized valence quark distributions in the 
proton as 


5u v (x, /i 2 ) = 40.3 x 2 ' 85 (1 — a;) 2 ' 15 , 


(2.2.4) 



and 


5d v (x,y 2 ) = —18.22a; 1 ' 41 (l — a;) 4 ' 0 . 


(2.2.5) 


Sg(x, y 2 ) = 0, 6qi(x, y 2 ) = Sq^x, y 2 ) = 0, (2.2.6) 

they are plotted in Fig.l. 

We predict the polarized glnon distribution at different Q 2 in Fig. 3. The results 
clearly show the accumulation of polarized gluons at small x. 

There are several databases of the polarized parton distributions, which are extracted 
by the global fitting DIS data. For example, we compare our results with the GRV 
distribution [22] in Fig. 4. The difference is obvious. It is not surprise that the polarized 
gluon distribution has large uncertainty, since the shape of the input gluon distribution 
is not constrained well enough by the DIS data alone. 

In order to understand the contribution of large gluon polarization, we draw the evolu¬ 
tion kernels P qg (z), P gg (z ), A P qg (z), AP gg (z) and yP gg ^ g (z),yP gg ^ q (z), yAP gg ^ g (z),yAP gg ^ q (z) 
in Figs. 5 and 6. One can find that 

(i) Pgg(z) > 0 and A P gg (z) > 0 imply that 5g(x, Q 2 ) is positive in our dynamic model; 

(ii) Since A P qg (z) < 0 at small z, we have dgi(x, Q 2 )/d\i\Q 2 ~ —Sg(x,Q 2 ) at small 
x, i.e., a large positive 5g at small x is expected to drive g\ towards large negative values; 

(iii) A Pgg^g < 0 and A P gg ^ q < 0 lead the net positive corrections of the gluon fusion 
to the polarized parton distributions since a negative sign in the shadowing terms of Eq. 
(2.1.3). To illustrate this effect, in Fig.7 we present xSg(x,Q 2 ) with and without the 
corrections of gluon recombination corrections at Q 2 = 1 and 5 GeV 2 . One can find that 
the effects of the gluon recombination in the polarized gluon distribution is positive. 

As we know that some approaches are planned to measure the gluon distributions. 

For example, the semi-inclusive deep inelastic scattering processes measure the 8g/g from 
hclicity asymmetry in photon-gluon fusion. The COMPASS collaboration [23] have used 
this method and found a rather small value for 5g/g = 0.024±0.080±0.057 at x = 0.09 and 
Q 2 = 3 GeV 2 . However we think that although the value of 5g/g is small, the polarized 
gluon contribution to the spin of the nucleon may be sizable since g itself is large at 
small x. In order to compare with the data, one needs to assume a suitable form for 
the unpolarized gluon distribution g(x,Q 2 ). Fortunately, both 5g(x } Q 2 ) and g(x,Q 2 ) are 
calculated within a same dynamics in this work and we avoid a larger uncertainty in the 
determination of g(x,Q 2 ). We compare our predicted 5g/g with the COMPASS data in 
Fig. 8. 

The other direct probing of 5g is offered by jet and 7r production in polarized proton- 
proton collisions available at BNL Relativistic Heavy Ion Collider (RHIC). A recent DSSV 
analysis [24] of high-statistics 2009 STAR [25] and PHENIX [26] data showed an evidence 
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of non-zero gluon hclicity in the proton. They found that the polarized gluon distribution 
in the proton is positive and away from zero in 0.05 < x < 0.2, although the presented 
data has very large uncertainty at small x region. Figure 9 presents the comparisons of 
our predicted gi(x, Q 2 ) at Q 2 = 10 GeV 2 with the DSSV bounds. Our results are beyond 
a up bound of the DSSV results, however, a sizable gluon polarization is still possible if 
taking the 90% confidence level (C.L.) interval. 

The NNPDF group has developed a new methodology [1] to extract polarized gluon 
distribution function. They used all essential available data and got an evidence of positive 
gluon polarization in the medium and small x region. This discovery is compatible with 
our results. Figure 10 shows the comparison of our predicted polarized gluon distribution 
with the NNPDF bounds. This example shows that a positive initial distribution of the 
polarized gluon at g 2 in the nucleon is impossible since it will obviously go beyond the 
up bound of the NNPDF analysis at Q 2 ~ 10 GeV 2 . 
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2.3 Discussions 


The total helicity of partons in a polarized proton are calculated by the first moments 

A f(Q 2 ) = f 1 dxdf(x,Q 2 ), (2.3.1) 

Jo 

Note that our predicted A q s (Q 2 ) for the sea quarks is positive since the negative contri¬ 
butions from the asymmetric strange quarks are neglected in this work. 

The contribution of quark polarization to the proton spin AE(Q 2 = 5 GeV 2 ) ~ 0.30 
in our estimation. It is interesting that this value is compatible with the world average 
values AE(Q 2 = 10GeO 2 ) = 0.31 ± 0.07 [27], and AE(Q 2 = 5 GeV 2 ) = 0.333 ± 0.011 ± 
0.025 ±0.028 [28], 

In Fig. 11 we plot the evolutions of AE(Q 2 ) and A g(Q 2 ) with increasing Q 2 . We 
find that A q v (Q 2 ) = 0.296, but A q s (Q 2 ) is slowly increasing from 0 at p 2 to 0.016 at 
Q 2 = lOOOGeV 2 due to the parton recombination corrections. On the other hand, the 
gluon helicity A g(Q 2 ) increases with In Q 2 beginning from zero at p 2 , and becomes large 
when Q 2 > 0.3 GeV 2 . 

The above mentioned AS and A g should be balanced by the orbital angular momenta 
of partons. For this sake, we write the nucleon helicity sum rule 

\ = 1a£(Q 2 ) + A g(Q 2 ) + y L‘(Q 2 ) + L;(Q 2 ), (2.3.2) 

where L z q denote the contributions of orbital angular momenta of quarks and gluons. In 
our model the sum rule at Q 2 = p 2 is 

I = IaE(P) + £LJ(P), (2.3.3) 

Z Z q 

where 

AS(p 2 ) = Au„(p 2 ) ± Ad y (p 2 ), (2.3.4) 

£L*(p 2 ) = L;(p 2 )±^(p 2 ). (2.3.5) 

<? 

On the other hand, the helicity sum rule of the polarized proton in its rest frame 
according to the constituent quark model is 

\ = \Y.^1 C ’ (2-3.6) 

where A q c denotes the quark polarization in the quark model. 

Comparing Eq. (2.3.6) with Eq. (2.3.3), we assume 
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1a«‘ = ±A «„(^)+Z4> 2 ), 

and 

+ (2.3.7) 

where the motions of partons are independent. 

Taking the SU(6) symmetry in the proton rest frame, we have 

Au c = 1 A <f = -i. (2.3.8) 

From Eqs. (2.2.5) and (2.2.6) we know 

A u v (fi 2 ) = 0.644, Ad v {fi 2 ) = -0.348. (2.3.9) 

Using Eq. (2.3.7) we obtain 

Z4(/i 2 ) = 0.345, L^(/i 2 ) = 0.007, (2.3.10) 

it implies that a polarized proton at scale /i 2 has two rotating n-qnarks, while the d-quark 
is located at the center of the proton since it has almost zero-orbital momentum. 

According to SU(6) symmetry, the constituent quark has zero angular momentum. 
However, according to Ref. [29,30,31,32] the orbital angular momentum L])(/i 2 ) in Eq. 
(2.3.7) may origin from the transverse distribution of the constituent quark in the rest 
frame due to the Melosh-Wigher rotation [33]. The Melosh-Wigher rotation is a pure 
kinematic effect in the frame transformation, we consider that this effect should keep the 
angular momentum conservation, for say, 

= b- s u + cL u , (2.3.11) 

where s u and L u are the spin and orbital angular momentum of a u-valence quark at / 1 2 in 
the light-cone frame of the polarized proton; 2/3 is from the SU (6)-distribution, the values 
of b and c depend on the wave function of valence quarks at fi 2 [29,30,31,32], Because 
the spin is an elemental physical quantity, it is always has |s£| = |s u | = 1/2, therefore, 
the orbital angular momentum L u changes only the direction of the spin from (it is 
also the polarized direction of the proton) to s u . Under these constraint conditions, once 
the values of b and c are determined by the wave function of valence quarks, the coupling 
angle between s u and L u in Eq. (2.3.11) can be fixed (see Fig. 12), and it leads to 
Au v (/i 2 ) < Au c . 
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We discuss the evolution of the sum rule (2.3.2) with Q 2 . The evolution equation for 
the quark and gluon orbital angular momenta at the leading order approximation was 
derived by Ji, Tang and Hoodbhoy in [34], it reads 


dE q L z q (t ) a s (t) 4 




a s (t) 


dt 


2tt 


dL Z g{t) 

dt 


Oisjt) 

2vr 


tc F J2L‘(t) 
6 <? 


Y l ‘M + 


a s (t) 

2n 


l-gC r S(t) 


where Cp = 4/3, rif is the number of active quark flavors, t 
ln(p 2 /A 2 ). The solutions are 


-Ag(t)], (2.3.12) 

ln(Q 2 /A 2 ) and t 0 = 


Eqw 2 ) 

q 

and 


-\mQ 2 )+ 


1 3 n f 

2 16 + 3 rif 


to q Z 


1 3n/ 

216 + 3n/ 


L^(g 2 ) = -Ag(Q 2 ) + 




16 


2 16 + 3n/ ’ 
(2.3.13) 


where /3 0 — 11 — 2n/3. Because of AE(g 2 ) = 0.296, L z g (fi 2 ) = 0, Ag(p 2 ) = 0, we can fixed 
^(A 2 ) = 1/2 - AS(p 2 )/2 = 0.352. 

Table. The contributions of various components to the proton spin at different Q 2 . 


Q 2 

0.064GeE 2 

1 GeV 2 

10GeE 2 

100GeE 2 

iAS 

0.148 

0.149 

0.151 

0.153 

E q L‘ 

0.352 

0.124 

0.096 

0.080 

A g 

0 

1.056 

1.993 

2.889 

L l 

0 

-0.829 

-1.74 

-2.622 

Total 

1 

2 

1 

2 

1 

2 

1 

2 


In Fig. 11 we add the curves of J2 q L Z (Q 2 ) and L g (Q 2 ) with rif = 3. These results 
describe a following novel spin-orbital picture of the proton in the light-cone frame: The 
proton is mainly constructed by one d-valence quark and two u-valence quarks, the d-quark 
is located at the center of the proton, and the two u-quarks rotate with J2 q L q (g 2 ) — 0.35 
around the d-quark at a bound state scale, with the Q 2 increasing the valence quarks 
radiate gluons and then sea quarks follow them. The former builds fast rotating glue 
cloud (see Table ), but their rotating direction is opposite to the u-quarks. Remind that 
the above mentioned possible orbital angular momentum of the partons in the polarized 
proton which can be checked in the experiments [35,36]. 
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3 Spin structure function g{ at small x 

3.1 A general consideration of the nucleon structure function 

at low Q 2 


In the research of the nucleon structure functions at the full kinematic region, an 
argued question is whether the parton distributions and their perturbative QCD evolution 
can (even partly) be applied to the low Q 2 range or the parton concept is suddenly invalid 
at a critical value of Q 2 < lGeV 2 ? 

Let us begin from the parton model for the spin-dependent distribution, which is writ¬ 
ten based on the Collins-Soper-Sterman (CSS) factorization schema [37] at the collinear 
approximation and in the twist-2 level, 

gi(x, Q 2 ) = / ( —Y / C q (x/y,Q 2 /fi F )Sq(y,fi F ), (3.1.1) 

J 0 y q 

which breaks up the spin structure function into two factors associated with perturbative 
short-distance functions C a and nonperturbative polarized parton distributions 5q at the 
factorization scale fi F . 

Taking the lowest order of C q 

C q (x/y, Q 2 /ii f ) = \e 2 q 8(x/y - 1 )S(Q - y F ) + 0{a s ) + 0(1/Q), (3.1.2) 

0(a s ) and 0(1/Q) are the QCD radiative corrections and higher twist contributions. 
Inserting it to Eq. (3.1.1), we obtain the relation between the spin structure functions 
and the polarized quark distributions 


Si(x,Q 2 ) = \Y.4Mx,Q 2 ) + iq(x,Q 2 ))+0(a,) + 0(llQ). (3.1.3) 

2 <7 

According to the renormalization group theory, 


dgi(x,Q 2 ) 
d In /ip 


(3.1.4) 


it gives the DGLAP equation 

Q 2 j 7 ^dq(x,Q 2 ) = I — P qq >(x/y,a s (Q 2 ))8q'(y,Q 2 ), (3.1.5) 

dQ Jo y ^ 

A P qq t denotes the splitting functions. If we consider only the leading order ( LO) approx¬ 
imation, we have 


g?GLA p ( x ,Q 2 ) = i^eJ[<5g(a;,Q 2 ) + 8q( x ,Q 2 )\, (3.1.6) 

z q 
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These results are available at Q 2 > a few GeV 2 . 

At lower Q 2 , the multi-parton correlations are important and the inclusive lepton- 
nucleon cross section is dominated by complicate higher twist terms. In fact, according to 
the operator product expansion (OPE), the spin structure function in the proton g\(x , Q 2 ) 
can be expressed as a series in 1 /Q 2 , 

9 i(x, Q 2 ) = gi T (x, Q 2 ) + g? T (x, Q 2 ). (3.1.7) 

The leading (twist-2) term corresponds to scattering from a single free parton, while higher 
twist terms correspond to multi-parton interactions. Only a little of higher twist can been 
calculated perturbatively in terms of quark and gluon degrees of freedom. For example, the 
contributions of parton recombination at initial (or finite) state to the DGLAP evolution 
equation have been calculated at leading order [10,11,12,13] and we denote this result as 
gDGLAP+ZRS ^ Q- 2 y However, we can neither perform nor interpret a partonic calculation 
of the higher twist effects containing the correlations between the initial and finite partons 
since they break the factorization schema. In a certain kinematic regime, some of such 
higher twist contributions to gi(x, Q 2 ) appear as observable hadronic phenomenon. In 
this case, we may chose a suitable phenomenological model, even do a parametrization to 
describe the corresponding higher twist effects. 

We try use the well known Vector Meson Dominance (VMD) model to mimic the 
above mentioned higher twist corrections. The reasons are as follows. The handbag 
diagram Fig.13a is a typical time ordered diagram describing Eq. (3.1.1), where the quark 
propagators connect with the probe and the target has only the forward components and 
these propagators can be broken as shown in Eq. (3.1.1) since they are on-mass-shell. 
The corresponding backward quark propagators construct the cat’s ear diagram Fig. 13b, 
which are neglected since these backward propagators are absorbed by the target in the 
collinear approximation [38,39]. However, the contributions of Fig. 13c can not been 
neglected at low Q 2 due to the corrections of quark-antiquark pair, which interacts with 
the target as a virtual vector meson if the transverse momentum k± ~ Q of quark pair 
is not large and confinement effects are essential. The contribution of Fig. 13c can not 
factorized as eq. (3.1.1). We use a phenomenological VMD model [40,41,42] to ’’isolate” 
this contribution from Fig. 13c. Traditionally, such VMD model was used to explain the 
structure function at low Q 2 region [43,44], We denote this contribution as g\ MD (x,Q 2 ). 

The more complicated corrections to g i at low Q 2 are from the higher order QCD 
effects O(aj-) and higher order recombination. In principle, we need to consider all these 
contributions, while it’s beyond our ability. Our motivation is that if one finds empirically 
that higher order corrections are deduced with a suitable scale down to low Q 2 (~ g 2 ), 
then one can extend our leading order analysis of structure function data to g 2 . If, the 
results are incompatible with the data, then the data can be used to extract the higher 
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order contributions. 

In consequence, at lower Q 2 we have 

9i(*, Q 2 ) ~ Pg? aLAF * ZRS (x, Q 2 ) + g FMD (x, Q 2 ) + g» T (x, Q 2 ), (3.1.8) 

where P is the probability of inelastic events via bare photon-parton interaction, the 
last term is the remaining higher twist corrections and we will neglect it at small x. 
Equation (3.1.8) implies that although the polarized partons share all nucleon’s spin, the 
higher twist effects mix with the contributions of partons in the measuring spin structure 
function at low Q 2 . We emphasize that g\ AID (x, Q 2 ) and gi T (x,Q 2 ) are irrelevant to 
the definition of the parton distributions because they violate the factorization schema, 
therefore, their contributions to gi(x,Q 2 ) will not change the discussions about spin in 
our previous section, which are the results of the polarized parton distributions in the 
proton. 
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3.2 Contributions of parton distributions and VMD part 


The contributions of the polarized parton distributions of the proton to the spin struc¬ 
ture functions at low Q 2 are 


2 ) = +Sq( X , g 2 )]. (3.2.1) 

2 q 

We assume that all parton distributions are freezed at scale Q 2 if Q 2 < p 2 . Based on this 
assumption we avoid the un-physical singularities at Q ~ A qcd- 

We present rr-dependence of gX GLAP+ZRS (x,Q 2 ) at several values of Q 2 in Fig. 14. 
One can find the dramatic change of the spin structure function at x < 10 ~ 3 from a flat 
form to dramatically decreasing. Considering Fig. 3, we conclude that the large gluon 
hclicity effect leads to this phenomenon. 

As we have mentioned that the contribution from the vector meson in virtual photon 
to g\ at Q 2 < 1 GeV 2 is necessary. According to the VMD model [43,44], 


xg\ MD 


(x,Q 2 ) 


1 rripQ 2 

8n-f 2 (Q 2 + ml) 2 




(3.2.2) 


where y p is the coupling constant of p vector meson and proton; We consider the contri¬ 
butions of p meson since 7 P <C 7 ^ <C 'yy, x is a variable defined as x = Q 2 /(s + Q 2 — m 2 ) 
rather than a momentum fraction of parton, s is the CMS energy square of the 7 p colli¬ 
sion. The cross-sections Aa pp (s) is the total cross section for the scattering of polarized 
meson with the nucleon, unfortunately, they are unknown. Usually, the following Regge 
theory [45] is used, 


A<j pp (s) ~ s x 1 ,ats—> 00 . (3.2.3) 

The extrapolation of g\ from the measured region down to x ~ 0 suggests us to assume 
that A = 1 — e and e ~ 0 is a small positive parameter due to the requirement of 
integrability of g± at x —> 0. In this work, we take e = 0. Thus, we have 


gX MD {x,Q‘ 2 ) ~ B- 


m 2 Q 2 


x X (1 — xy, 


(3.2.4) 


0 Q 2 + rn 2 p ) 2 

where B = 0.03 and the factor (1 — x) 7 is due to the spectator counting rules at high x 
[46], and it restricts the application of the VMD model in small x range. 

We read m A /(Q 2 + m 2 ) 2 in Eq. (3.2.4) as the probability of the VMD event, therefore, 


P = 1 - 


rn; 


0 Q 2 + ml) 2 ' 


(3.2.5) 


in Eq.(3.1.8). 
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3.3 Predictions for spin structure function g\ at small x 


What is the asymptotic behavior of This is a broadly discussed subject. We plot 
gKx^Q 2 ) = Pg? GLAP+ZRS (x,Q 2 ) + gY MD (x,Q 2 ) with different values of Q 2 in Fig. 15. 
There are two different asymptotic behaviors of gj 3 at small x: the VMD behavior ~ x^ 1 
at Q 2 < 1 GeV 2 and the large gluon helicity effect at Q 2 > ?>GeV 2 . Besides, g{ presents 
the twist form of the two asymptomatic behaviors above, which is the mixing result of 
the nonperturbative and perturbative dynamics. 

We compare our predicted g T { at x > 10” 3 with the data [47] in Fig. 16. These data 
on 2010 are more precise than the previous data. Note that the values of Q 2 of every 
measured point are different and they are taken from Table I of [47]. The theoretical curve 
is a smooth connection among these points. This figure shows that the pQCD evolution 
almost control the behavior of g\ at x > 10 -3 . 

On the other hand, the combination of nonperturbative and perturbative dynamics at 
x < 10 -3 leads to a dramatic change of g{ around Q 2 = 1 ~ 3 GeV 2 . Unfortunately, there 
are only several data with large uncertainty about gj 3 in this range. In Figs. 17 and 18 we 
collect the HERA early data [48,49] at Q 2 = 1, lOGeU 2 which are un-generally used and 
compare them with our predicted g[. Figure 19 shows some of these data (trigon) and the 
comparisons with our results (dark points). Figure 20 is the Q 2 -dependence of g\ with 
fixed x , the data are taken from [50]. One can find that our predicted g j 3 are compatible 
with these data, although more precise measurements are necessary. 

Finally, we compare our results with the new COMAPSS (primary) data [14,15,16,17,18] 
at Q 2 < 1 GeV 2 , which show that gj 3 presents a flat asymptomatic form at x < 10 -3 . This 
seems to contradict with the predicted strong rise of g\ at Q 2 < 1 GeV 2 in Fig.3. However, 
in the COMPASS fixed target experiments there is a strong correlation between x and Q 2 , 
which makes it possible that low x measurements are along with low Q 2 . In Fig. 21, we 
take the average values of Q 2 for each probing values of x (see Fig.l in Ref.[15-18]). The 
results are acceptable. Obviously, the measurements at different x with different values 
of Q 2 in the fixed target experiments mix two different asymptomatic behaviors of g\. 

We predict the stronger Q 2 - and ^-dependence of g{ at 0.01 < Q 2 < 3 GeV 2 and x < 0.1 
due to the mixture of nonperturbative vector meson interactions and the QCD evolution 
of the parton distributions in Fig. 15. For testing this prediction, the measurements of g\ 
with fixed x or Q 2 at low Q 2 are necessary. The planning Electron-Ion Collider (EIC), 
for example, eRHIC [51] and EIC@HIAF [52] can probe a broad low Q 2 < lGeU 2 -range, 
where we can check the predicted behavior of g{ at fixed x or Q 2 . 
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3.4 Discussions 


In general consideration, both the logs of \ jx and Q 2 are equally important at small 
x and low Q 2 , and one should sum the double logarithmic (DL) terms (ct s ln 2 (l/a;)) n , 
which predict the singular behavior g\ ~ x~ x (A > 0). It means that the BFKL equa¬ 
tion [53,54,55,56,57,58] and its nonlinear corrections- the Balitsky-kovchegov equation 
[59,60,61] and the JIMWLK equation [62,63,64,65,66,67] should combine with the DGLAP 
equation. However, the translation between the BFKL equation and the DGLAP equation 
is a complicated technic. One of such method is the Ciafaloni-Catani-Fiorani-Marchesini 
(CCFM) equation [68,69,70,71], which is derived based on the two-scale unintegrated 
gluon distribution. The solution of the CCFM equation is much more complicated and 
has only proven to be practical with Monte Carlo generators. To avoid this difficulty, 
some special methods are proposed [72,73]. For example, the double logarithmic terms 
are taken into account via a suitable kernel of the evolution equations in the infrared 
evolution equations, which was first suggested by Lipatov [74,75], or alternatively taking 
a singular initial parton distributions at x < ICC 2 , one can also mimic the results of the 
DL-resummation. 

In this work, the behavior of gi at the same range is obtained through a long evolu¬ 
tion of the DGLAP equation with the parton recombination corrections. We ford that it 
is different from the predictions of the DL-resummation, the asymptomatic behavior of 
the polarized quark distributions at x —> 0 is controlled by A P qg in the DGLAP equa¬ 
tion, rather than the ln A (l/a;)-corrections to the DGLAP-kerncl. Thus, the difficult DL 
resummation can be replaced by the fits of the initial quark distributions 8q v (x,/j> 2 ) in 
the DGLAP equation if the evolution distance is long enough. This conclusion was also 
obtained in the unpolarized structure functions [76]. 
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4 Origins of the generalized Gerasimov-Drell-Hearn 

sum rule 

4.1 Spin structure functions in the full Q 2 range 


The Gerasimov-Drell-Hearn (GDH) sum rule reads 

o M 2 k 2 

m = lim r !(g 2 ) = -- ~ -0.8, (4.1.1) 

Q 2 ->0 4 

where k is the anomalous magnetic moment of the nucleon. On the other hand, the 
Bjorken sum rule [77] says 


lira [r?(Q 2 ) - r;(Q 2 )] = 1 

Q 2 ->oo 6 


9 A 


9v 


(4.1.2) 


this ratio is accurately known as[78]: qa!9v — —1.2695 ± 0.0029. 

The connection of the two sum rules by means of the generalized GDH sum rule is(For 
an overview, see Ref. [79] for example) 


9 M 2 

h(Q 2 ) = -grUQ 2 ), (4.1.3) 

which allows us to study the transition between the perturbative partonic structure and 
nonperturbative hadronic picture of nucleon in lepton-nucleon scattering processes. The 
data show that this sum rule at low Q 2 < IGeV 2 changes dramatically and exceeds the 
variation bound at higher Q 2 , which has been parameterized (but not explanation) in 
[80,81,82,83]. The explanation of the generalized GDH sum rule is an active subject. For 
example, the phenomenological constituent quark model [84,85], the VMD model [86,87], 
the resonance contributions [88], the chiral perturbation theory (yPT) [89,90,91] are used 
to understand the generalized GDH sum rule. 

The first moment of g\ is 


r?(Q 2 ) = y° glap+zrs (q 2 ) + y\ md {Q 2 ) + rf T (Q 2 ) 

From Eqs. (3.1.8) and (3.2.5) we obtain 


■^DGLAP+ZRS ( r> 2 


(Q 2 ) ~ 0.123 1 - 


ml 


(4.1.4) 


(4.1.5) 


0 Q 2 + ml) 2 J ' 

The dashed curve in Fig. 22 is our predicted Y PGLAP+ZRS (Q 2 ). On the other hand, 


we have 


T MD (Q 2 ) = I 1 dxgX MD {x,Q 2 ) ~ 0.055 
Jo 


m 2 v Q 2 


0 Q 2 + ml) 2 ' 


(4.1.6) 
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Comparing the solid curve with T^Q 2 ) data [92,93,94,95,96,97,98] in Fig. 22, one can 
expect that the remaining higher twist corrections play a significant role at low Q 2 
to the general GDH sum rule. We will discuss them in detail next section. 
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4.2 Higher twist contributions to the GDH sum rule 


According to the OPE, the appearance of scaling violations at low Q 2 is related to the 
higher twist corrections to moments of structure functions. Higher twists are expressed 
as matrix elements of operators involving nonperturbative interactions between quarks 
and gluons. The study of higher twist corrections gives us a direct insight into the nature 
of long-range quark-gluon correlations. The higher twist corrections to (j\ have several 
representations. In this work, we will try to expose the remaining power suppression 
corrections to Tf. For this sake, we make Y^ T {Q 2 ) (i.e., the data points in Fig. 22)- 
^dglap+zrs^q 2 ) + Figure 23 shows such a result at Q 2 > 0.2 GeV 2 , which 

has been smoothed with minimum p/D.o.f.. 

To expose the possible physical information of the curve in Fig. 23, according to QCD 
operator product l/g 2 -expansion, 


-pHT/^2\ _ y- h-2 i{Q 2 ) 

1 1 Q2i~2 ’ 

we take first three approximations 


rf T<4| W 2 ) = 

r HT(4+6) ,_ h6 (Q 2 ) + ^a{Q 2 )Q 2 

1 1 \ye ) ~ 7u 


(4.2.1) 


(4.2.2) 

(4.2.3) 




HT{ 4+6+8)/ n 2N _ 


m = 


hs (Q 2 ) + ^e(Q 2 )Q 2 + i^a{Q 2 )Q a 

eft 


(4.2.4) 


Then we plot the curves Q 2 rf T(4) (Q 2 ), Q 4 rf T(4+6) (Q 2 ) and g 6 rf T(4+6+8) (g 2 ) in Fig. 
24. There are following interesting properties of these results: 

(i) g 6 rf T ( 4 + 6 + 8 )(g 2 ) —). 0, if g 2 —> 0. This implies that vanishes if it is independent 
of g 2 . Therefore, Tf 


HT( 4+6)/ n 2 


(g 2 ) is an appropriate approximation. 

(ii) Three curves in Fig. 24 cross at a same point Q 2 ~ 1 GeV 2 . Particularly, the 
intercept /jq of the line suddenly changes its value from -0.037 at Q 2 > 1 GeV 2 to 0.006 
at Q 2 < 1 GeV 2 . This result exposes that the correlation among partons in the proton 
has an obvious change near Q ~ 1 GeV. 

(iii) We use 


rf T(4+6) (g 2 ) = 


Hi 


+ 


h 6 


at Q 2 < 0.3 GeV 2 


(4.2.5) 


g 2 + e 2 (g 2 + e 2 ) 2 

to fit the data at Q 2 < 0.3GeV 2 , where we add a parameter e to remove the unnatural 
singularity at Q 2 = 0. The value of e is sensitive to /(0). We find that ji 4 = —0.13GeP 2 , 
= 0.0528GeF 4 and e 2 = 0.422GeF 2 . 
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In summary, 


T?(Q 2 ) = 0.123 1 


ml 


(<Q 2 + ml ) 2 / 

where the HT contributions are 


+ 0.055 


m 2 v Q 2 


0 Q 2 + m 2 v ) 2 


+ rf T (g 2 


(4.2.6) 


rf T (Q 2 ) = 


0.004M 2 0.037M 4 

Q-> n 3 

0.048M 2 __ 

Q 4 
■ 4 - 


O 2 
0.13M 2 


•i 9 i n i no 71 r 9 


— at Q 2 > 1 GeV 2 
+ °'°^ M4 at 0.3 < Q 2 < 1 GeV 2 

0.0528M 4 . ri2 ^ 


~v~ ,-Q 


(4.2.7) 


where M 2 = lGeV" 2 . We present the comparison of our Tf(g 2 ) with the data [92,93,94,95,96,97,98] 
in Fig. 25. The corresponding /f(Q 2 ) is presented in Fig. 26. 
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4.3 Discussions 


The parton-hadron duality was first noted by Bloom and Gilman [99,100] in deep 
inelastic scattering (DIS) and has been confirmed by many measurements. At low ener¬ 
gies (or intermediate Bjorken variable x and low Q 2 ) DIS reactions are characterized by 
excitation of nucleon resonances; while at high virtuality such processes have a partonic 
description. The smooth high-energy scaling curve essentially reproduces the average of 
the resonance peaks seen at low energies. Burkert and Ioffe [81] indicated that the con¬ 
tribution of the isobar A(1232) electro-production at small Q 2 can describe the general 
GDH sum rule, and they gave 

= -0.056 ~ -0.063, at Q 2 = 0.3 ~ 0.8 GeV 2 (4.3.1) 

M z 

= 0.010 ~ 0.011, at Q 2 = 0.3 ~ 0.8 GeV 2 , (4.3.2) 

which are compatible with our prediction Eq. (4.2.7). 

Our results in Fig.24 indicate that the negative twist-6 and twist-4 effects dominate 
the suppression of Tf(Q 2 ) at Q 2 > 1 GeV 2 and Q 2 < 1 GeV 2 . Particularly, the slope of p 4 
of the lines, which cuts Q 4 rf T( ' 4+6 \Q 2 ) suddenly changes its sign at Q 2 < 1 GeV 2 . This 
result exposes that the correlation among partons in the proton become stronger at scale 
~ 0.2 f m. We noted that Petronziol, Simula and Ricco [101] reported that the inelastic 
proton data obtained at Jefferson Lab exhibit a possible extended objects with size of 
~ 0.2 — 0.3 fm inside the proton. 
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5 Summary 


In this work we consider that the nucleon is consisted of quarks and gluons (partons) 
via QCD interactions even at low Q 2 . A general worry is that the correlations among 
partons may break the definition of parton distributions and their evolution rules. As 
a model, we treat these high twist effects as two parts: (i) The leading recombination 
among initial partons, which modifies the DGLAP equation but keeps the momentum 
conservation in Eqs. (2.2.1)-(2.2.3) and the nucleon hclicity sum rule Eq. (2.3.2); (ii) 
The phenomenological VMD model and the parameterized higher twist effects, which 
contribute to the measured structure functions of the nucleon but they are irrelevant 
to the parton distributions in the nucleon. In this framework, we discuss the electron 
scattering off a nucleon at high energy in a special (Bjorken) infinite momentum frame, 
where the virtual photon presents two components: bare photon 7* and vector meson 
V*with J PG = 1 . In the former case, 7* couples either with an on-mass-shell quark 

and contributes p PGLAP + ZRS ; where we take a leading order approximation and all higher 
order corrections are absorbed into the free parameters, or with an off-mass-shell quark, 
which gives F PT . In the later case, the VMD model describes the nonperturvative multi- 
parton interactions between V* and nucleon. Thus, we present a compact theoretical 
model about the nucleon spin structure. 

(i) We find that the gluon contribution to the spin of proton is much larger than the 
predictions of most other theories. This result is compatible with the recent NNPDF 
analysis and suggests a significant orbital angular momentum of gluons to balance the 
contribution of gluon spin. In concretely, the total proton spin at a bound state scale g 2 
is composed by ~ 30% quark spin and ~ 70% orbital angular momentum of the quarks, 
where two u-valence quarks are rotating around a d-valence quark. With increasing Q 2 , 
the omitted gluons accumulate a larger positive helicity, which is mainly balanced by their 
orbital momentum. Therefore, there are two rotating groups in a polarized proton at Q 2 \ 
a slower quark group and a faster gluon. 

(ii) We use the DGLAP equation with the parton recombination corrections and the 
nonperturbative VMD model to predict the spin structure functions g\ of the proton. We 
first present a complete picture for the translation of g\ from low Q 2 (~ 0) to high Q 2 at 
small x. We find that the contribution of the large gluon hclicity dominates g\ at x > 10 -3 , 
but the mixture with nonperturbative component complicates the asymptomatic behavior 
of g\ at x < 10 -3 . The results are compatible with the data including the early HERA 
estimations and COMPASS new results. The predicted strong Q 2 - and x-dependence of 
g\ at 0.01 < Q 2 < 3 GeV 2 and x < 0.1 due to the mixture of nonperturbative vector 
meson interactions and the QCD evolution of the parton distributions can be checked on 
the next Electron-Ion Collider (EIC). 

(iii) We discuss the contributions of parton distributions and VMD component to the 
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lowest moment of the spin-dependent proton structure function. After removing the above 
two contributions from the existing experimental data for rf(Q 2 ), the higher twist power 
corrections present their interesting characters: parton correlations at Q 2 rsj 1 GeV 2 show 
a bend point, where the twist-4 and twist-6 effects dominate the suppression of Y\(Q 2 ) 
at Q 2 < 1 GeV 2 and Q 2 > 1 GeV 2 , respectively. The results suggest a possible extended 
objects with size 0.2 — 0.3 fm inside the proton. Within the analytic of these results, we 
are able to achieve a rather good description of the data at all Q 2 region using a simple 
parameterized form of Ti(Q 2 ). 
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Appendix: 

From Ref. [13], we have LO polarized gluon recombination functions 

Pg+g+^g + = 7 -- -—-7T— - 2 x x 3 x + ay V + X 2 4 ~ 2 X 2 X + X 2 X 2 

4 XX 2 Z (ay + X 2 ) X\ z 

-\-X\X 2 — X\X 2 X — X\X 2 X + XiX 2 X 2 ) (A.l) 

9 (ay + X 2 — x) . 4 f /? 3 2 q 3 2 72 3 

Pg+g+^g- = 7-—-73—- (6 ay aya; + 6 Xi x 2 x -3 ay x 2 x - 7 X 1 X 2 X 

4 xx 2 z (X\ + X 2 ) X\ z 

+11 x 2 x 2 x 2 + 6 x 2 x\x + 6 x 2 x 2 x — 3 x 2 3 x\x 2 — 7 ay 2 aya; 3 + 2 ayayar 
+ay l 1 + x 2 + 2 ay 5 ay + 2 ay 4 ay 2 — ay 4 a; 2 — 2 ay 3 a; 3 + 2 ay 2 a; 4 + 2 ay 5 ay 
+2 ay 4 ay 2 — ay 4 a; 2 — 2 ay 3 a; 3 + 2 x 2 x 4 + 2 ay 3 ay 3 ) (A.2) 

Pg+g-~> 9 + = t — ^+7 + X2 —r=r—7(141 ay 7 a; 2 a ;2 + 42 ayVay 2 - 100ay 8 a;ay 
4« 2 2 (xi+a:2) ii 2 

+19 xi 5 x 4 x 2 — 5 ay 8 aya; + 39 x 2 x 4 x 2 + 137 ay 3 ay 6 at — 86 ay 6 a; 3 ay 

+26 x 2 9 xi + 5 ay 10 + 5 x 2 10 + 155 ay 4 a; 2 ay 4 — 40 x 2 x 2 x 3 — 212 x\xx 2 

—124ay 3 a; 3 a;2 4 — 196 ay 4 a? 3 ay 3 — 79 x 2 x 2 x 2 + 128ay 4 ay 5 a; + 44ay 7 ay 2 a; 

—47 ay 5 a?ay 4 — 177 x\x 3 x 2 — 209 x\ xx 2 + 40 ay 3 a; 4 ay 3 — 33 ay 3 ay 5 a; 2 

+ay 6 aya; 3 + 319 x\x 2 x 2 + 291 x\x 2 x 2 — 35 ay 7 aya; 2 + 13 ay 5 a; 4 ay 

+64 ay 7 ay 3 + 57ay 8 ay 2 + 26 ay 9 ay + 57ay 8 ay 2 + 64 ay 7 ay 3 + 34 ay 6 ay 4 

+12 ay 5 ay 5 + 34 ay 4 ay 6 — 4 x 2 x + 2 x 2 6 x A + 2 x 2 x 3 — 5 x 2 s x 2 + 30 ay 8 a; 2 

—20 ay 7 a; 3 + 5 ay 6 a; 4 — 20 ay 9 a;) (A.3) 

9 (ay x 2 — x) , 7 2 ,07 4 4 2 78 , n 5 4 

Pq+q-^q. = - -77-77—+-31 ay 'x x 2 + 27 Xi x x 2 - 7ay° xx 2 +9x7x x 2 

A xx 2 z (xi + x 2 yxi z 

— 104ay 8 aya; + 54ay 2 a; 4 a;2 4 — 206 ay 3 ay 6 a; + 3ay 6 a; 3 a;2 + 26ay 9 ay + 5 ay 10 
+5 ay 10 + 115 ay 4 a; 2 ay 4 — 211 x 2 x 2 x 3 + 125 x\xx 2 — 192 x 3 x 3 x 2 4 
—92 ay 4 a; 3 ay 3 + 319 x 2 x 2 x 2 — 25 ay 4 ay 5 a; — 217 x 2 x 2 x + 136 ay 5 a;ay 4 
—24ay 5 a; 3 ay 2 + 34ay 7 a;ay 2 + 36 ay 3 a; 4 ay 3 + 307 x 3 x 2 x 2 — 106 ay 6 aya; 3 
—41 ay 5 a? 2 ay 3 — 67 x\x 2 x 2 + 157 x 2 x\x 2 + 27ay 5 x 4 ay + 64 ay 7 ay 3 
+57 ay 8 ay 2 + 26ay 9 ay + 57ay 8 ay 2 + 64 ay'ay 3 + 34ay 6 ay 4 + 12ay 5 ay 5 
+34 ay 4 ay 6 — 20 x 2 x + 5 x 2 x A — 20 x 2 x 3 + 30 x 2 x 2 — 5 ay 8 a; 2 + 2 x \' x 3 
+2 a;i h a; 4 — 4 ay 9 a;) (A.4) 

Setting x\ = x 2 = y , one can find 

^Pgg^-g = [Pg+g+-*-g+ ~ Pg+g+^-g- ~ Pg+g--*g+ + Pg+g--*g-] 
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Similarly, 


p 

1 g+9+^q+ 


p 

r g+g+-+q- 


p 

- r 9+9--a?+ 


P 

1 9+9- -a?- 


Setting x\ 


= [P t 


— p 

9 + 9 +^ 9 + 1 g+ 9+^9 


+ Pa 


Pn 


g+g-^ 9 + ~ r g+g-^> 9 - 

27 (2 y — x)(—20y 3 + 12y 2 x — x 3 ) 

64 v 5 


■A / 11 9 (4xi 4 + 7 x- 3 x 2 -8x 1 3 x + 2x 2 x 2 -6 x 1 2 xx 2 

12 (#1 + x 2 ) 6 x 2 z xp 

+4 x 2 xi 2 + 4 x 2 x 2 2 — xix 2 3 + 2 x\xx 2 2 — x\x 2 x 2 ) (A.5) 

~ ^ 1 + — j^t(4x 2 x 1 2 + 4xx 2 x 2 2 + 8xxx 2 3 - 8xxxx 2 2 + 4x 2 

12 [xi + x 2 ) 6 x 2 2 xi 2 

—8x 2 3 x + 4r 2 x 2 2 — xix 2 x 2 ) (A.6) 

-*r } X 1 + — p-^(4x! 6 x 2 + 4 x! 4 x 2 4 - 8 x\ x + 24x\'x 2 2 + 16 xi 3 x 2 3 

12 (x 1 + x 2 y x 2 X\ 

+ 16 :ri 7 :?’o -I- 4 :ri 8 + 24 r 1 




—8 Xi 4 a;a: 2 3 + 24 xi 5 x 2 x 2 + 33 xi 4 x 2 x 2 2 + 14 rci 3 :^ 4 — 40 xi & xx 2 

—53 x\xx 2 — 8 xx 2 x 2 — 9 x\xx 2 ) (A.7) 

1 {xi + x 2 — (24 X2 6 Xi 2 + iq x J Xl - 8 x 2 x + 16a^i 3 + 4x 2 8 
12 (24 + £ 2 ) X 2 

+4 rc/A 2 + 4 rc/a^ 4 + 6 x 2 x 2 x\ + 4 a; 2 ^ 6 + 8 x 3 x 2 x 2 + 15 xi 2 a; 2 a: 2 4 
—13 x^xx 2 + 24 x 1 5 x 2 x 2 + 51 rri 4 :*; 2 ^ 2 — 35 + x\xx 2 

—35xa;2 5 a;i 2 — 22x 1 xx 2 6 ) (A.8) 


= x 2 = y, we have 


A P t 


gg^-q ~ [Pg+g+->q+ 


Pn 


g+g+^q- 

= [Pg+9+^-q+ ~ Pg+g+^rq- A Pg+g-^q+ ~ 1 g+g-^q 
1 (2 y-x) 2 (-y + x) 


Pg+g-^q+ + Pg+g-^q- 
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Figure 1: Input valence quark distributions at g 2 = 0.0 QAGeV 2 for polarized and un¬ 
polarized densities. 




Figure 2: Fitting the data [21] for g\ at Q 2 = 5 GeV 2 using the input distributions Eqs. 
( 2 . 2 . 1 )-( 2 . 2 . 6 ). 
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Figure 3: Predicted polarized gluon distribution x5g(x, Q 2 ) in the proton at the different 
Q 2 -scales, which show the accumulation of radiative polarized glnons at small x in the 
evolution. 






Figure 4: Comparisons of our predicted polarized LO parton distributions at Q 2 = 5 GeV 2 
with the GRV distributions [22], 
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Figure 5: The splitting functions. 
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Figure 6: The recombination functions. 
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Figure 7: Predicted polarized gluon distributions xSg in the nucleon at Q 2 = 1 and 5 GeV 2 
with gluon recombination corrections (solid curves) and without gluon recombination 
corrections (dashed curves). 



Figure 8: Comparison of dynamically predicted 5g/g with the COMPASS data [23] at 
Q 2 = 1 and 5 GeV 2 . 
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Figure 9: Our predicted x5g at Q 2 = 10 GeV 2 (solid curve). Broken curve is the result 
by DSSV using RHIC measurements in [24]; dotted curves are the fits within the 90% 
confidence level (C.L.) interval 



Figure 10: Comparison of dynamically predicted polarized gluon distribution with the 
NNPDF bounds [1]. 
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Figure 11: Contributions of spin and orbital motion of the partons to the proton spin and 
their evolutions with Q 2 . 



Figure 12: A schematic diagram of the proton spin crisis: Orbital angular momentum 
L u of a valence u-quark at a bound state scale /j, 2 impels the direction of the u-quark spin 
(s u ) to deviate the polarized direction of the proton (s£) and gives AS < 1. 
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Figure 13: The time ordered decomposing of DIS diagrams, (a) The struck quarks are 
on-mass-shell since they have only forward component, (b) A ’’cat ear” diagram, which 
vanishes in the collinear factorization schema, (c) The ’’cat ear” diagram with higher 
order QCD corrections, which are non-vanished at low Q 2 , but can be isolated using a 
naive VMD model. 


i 



Figure 14: Perturbative g P)GLAP + ZRS . All partons are evolved from three valence quarks 
at jj 2 = 0.0 64GeV 2 . 
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Figure 15: g\ evolutions at different values of Q 2 in (a) large and (b) small scales. 



Figure 16: Predicted g{ at x > 10 3 and comparisons with the COMPASS data [47]. Note 
that the values of Q 2 (x ) of each measured point are different (see Table I of Ref. [29]). 
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Figure 17: Predicted g\ at Q 2 = 1 GeV 2 and the comparison with the HERA data [48]. 



Figure 18: Predicted g\ at Q 2 = 10 GeV 2 and the comparison with the HERA ’’data”, 
which are based on the NLO QCD predictions with the statistical errors expected at 
HERA [49], 
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Figure 19: Predicted g[ at Q 2 = 1.8GeR 2 , 5.6GeV 2 and 16.5 GeV 2 at x < 10 3 (circles) 
and the comparison with the HERA data (triangles). 



Figure 20: Predicted (^-dependence of g\ with Exed valiiesof x. The data are taken from 

[50]. 
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Figure 21: Predicted g\ as a function of x with different measured Q 2 (x ) (solid curve). 
Note that the low values of x connect with the low values of Q 2 (x). The data are taken 
from COMPASS primary results with two different beam energies [14,15,16,17,18]. 


ft: 



Figure 22: Contribution of quark helicity y pglap+zrs (Q 2 ) (dashed curve) and combining 
VMD contribution y pglap+zrs (Q 2 ) + T] MD (Q 2 ) (solid curve). The data are taken from 
Hermes experiment at DESY [92,93], the E143 experiment at SLAC [94] and the EGla 
experiment using the CLAS detector at JLab [95,96,97,98]. 
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Figure 23: Contribution of higher twist T^ T (Q 2 ) (smoothed curve) is taken from data- 

[r DGLAP+ZRS ( q 2 ) + T VMD(Q 



Figure 24: Three different analysis of the higher twist contributions. 
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Figure 25: The Q 2 dependence of T^Q 2 ) calculated by Eqs. (4.2.6) and (4.2.7). The data 
are taken from [92,93,94,95,96,97,98]. 



Figure 26: The Q 2 dependence of Jf(Q 2 ) calculated by Eq. (4.1.3). The data are taken 
from [92,93,94,95,96,97,98], 
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